ANISOTROPIC THERMO-ELASTICITY IN 2D 
PART I: A UNIFIED TREATMENT 
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Abstract. In this note we develop tools and techniques for the treatment of anisotropic 
thermo-elasticity in two space dimensions. We use a diagonalisation technique to obtain 
properties of the characteristic roots of the full symbol of the system in order to prove L p —L q 
decay rates for its solutions. 
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1. The problem under consideration 

Systems of thermo-elasticity are hyperbolic-parabolic or hyperbolic-hyperbolic coupled sys- 
tems (type-1, type-2 or type-3 models) describing the elastic and thermal behaviour of elastic, 
heat-conducting media. The classical type-1 model of thermo-elasticity is based on Fourier's law, 
which means, that the heat flux is proportional to the gradient of the temperature. The present 
paper is devoted to the study of type-1 systems for homogeneous but anisotropic media in R 2 . 
There arc different results in the literature for certain anisotropic media (cubic in [3]; rhombic 
in [3]). Our goal is to present an approach, which allows to consider (an)isotropic models in R 2 
from a unified point of view. 

We consider the type-1 system of thermo-elasticity 

U a + A(D)U + 7 V6> = 0, (1.1a) 

6 t - nA6 + jV T U t = 0. (1.1b) 

Here A(D) denotes the elastic operator, which is assumed to be a homogeneous second order 
2x2 matrix of (pseudo) differential operators and models the elastic properties of the underlying 
medium. Furthermore, k describes the conduction of heat and 7 the thermo-elastic coupling of 
the system. We assume n > and 7^0. We solve the Cauchy problem for system (jl.ip with 
initial data 

U(0, ■) = U 1 , U t (0,-) = U 2 , 9(0,-) =0a, (1.2) 
for simplicity we assume Ui, U 2 € S(R 2 ,R 2 ) and a e <S(R 2 ). We denote by A(£) the symbol of 
the elastic operator and we set 77 = £/|£| G S 1 . Then some basic examples for our approach are 
given as follows. The material constants are always specified in such a way that the matrix A(j]) 
becomes positive. 

Example 1.1. Cubic media in 2D are modelled by 

yl) \ (A + fi) mm (r - fi)^ +fij v ' 

with constants r, [i > 0, — 2fi— r < A < r. This case was treated e.g. in [3]. For the corresponding 
elastic system see [T2"] . 

Example 1.2. Rhombic media in 2D are modelled by 

m = ( ( ?r + M , (A + ^V 2 1 (1.4) 
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with constants t\,T2, fx > and — 2[i — ^Jt\T2 < A < Jt\t^. For this case we refer also to @]. 

Example 1.3. Although it is not the main point of this note, we can consider isotropic media, 
where 

A{rj) = fxl + (A + fi)rj ® 77 

_ /(A + m)»7i + /* (A + ^)'7i^2 \ ^ 5 <, 

with Lame constants /i > and A + /1 > 0. 

We will present a unified treatment of these cases of (in general) anisotropic thermo-elasticity. 
For this we assume that the homogeneous symbol A = A(£) = \£,\ 2 A(n), r\ = £/|£|, is given as a 
function 

A : S 1 -» C 2x2 (1.6) 

subject to the conditions 

(Al): A is real-analytic in rj £ S 1 , 

(A2): A(ry) is self-adjoint and positive for all 77 € S 1 . 
For some results we require that 

(A3): A{r\) has two distinct eigenvalues, #specA(^) = 2. 
Under assumption (A3) the direction 77 <E S 1 is called (elastically) non- degenerate. In this case we 
know that the elasticity equation Utt +A(D)U — is strictly hyperbolic and can be diagonalised 
smoothly using a corresponding system of normalised eigenvectors Tj (77) to the eigenvalues *Cj (77) 
of A{rj). 

If (A3) is violated we will call the corresponding directions 77 £ S 1 degenerate. For these 
directions we can use the one-dimensionality of S 1 in connection with the analytic perturbation 
theory of self-adjoint matrices (cf. [6]). So we can always find locally smooth eigenvalues >Cj{rf) 
and corresponding locally smooth normalised eigenvectors rj(n) of A{rj). For the following we 
assume for simplicity that these functions extend to global smooth functions on S 1 . 

This classification of directions is not sufficient for a precise study of L p -L q decay estimates 
for solutions to the thermo-elastic system. It turns out that different microlocal directions 
77 = £/|£| £ S 1 from the phase space have different influence on decay estimates. But how to 
distinguish these directions and how to understand their influence? In general, this can be done 
by a refined diagonalisation procedure applied to a corresponding first order system (first order 
with respect to time). Applying a partial Fourier transform and chosing a suitable energy (of 
minimal dimension) this system reads as Y) t V — B(£)V. The properties of the matrix B(£) are 
essential for our understanding: 

• the notions of hyperbolic and parabolic directions depend on the behaviour of the eigen- 
values of B{£) (see ([2~3]) . 1(23)1 ■ Definition (J); 

• the matrix B{^) contains spectral data of A(£) together with certain coupling junctions 
(see (|2.6|) ) between different components of the energy. The behaviour of these coupling 
functions close to hyperbolic directions has an essential influence on decay rates (see 
Theorems O [33] and E3J . 

It turns out that we have to exclude some exceptional values of the coupling constant 7 by 
assuming that (see Definition [1] for the notion of hyperbolic directions) 

(A4): 7 2 7^ 2xj (fj) — tr A(fj) for all hyperbolic directions fj with respect to Xj . 
Basically this implies the non-degeneracy of the 1-homogeneous part of B(£). In the following 
we will call a hyperbolic direction violating (A4) a ^-degenerate direction. Assumption (A4) is 
used for the treatment of small hyperbolic frequencies and plays there a similar role like (A3) 
for large frequencies. 
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In Section [2] we will give the transformation of the thermo-elastic system (jl.ll) to a system 
of first order and the diagonalisation procedure in detail. The proposed procedure generalises 
those from [H] , [H] , [§] , [I7J . The obtained results are used to represent solutions of the original 
system as Fourier integrals with complex phases. Based on these representations we give micro- 
localised decay estimates for solutions in Section [31 They and their method of proof depend 
on 

• the classification of directions (to be hyperbolic or parabolic) ; 

• the order of contact between Fresnel curves (coming from the elastic part) and their 
tangents for hyperbolic directions; 

• the vanishing order of the coupling functions in hyperbolic directions. 

Let us formulate some of the results. The first one follows from Theorem [O] and Corollarv l3.2l 

Result. Under assumptions (Al) to (A4) and if the coupling functions vanish to first order in 
hyperbolic directions the solutions U(t,x) and 9(t,x) to |HTj) satisfy the L p -L q estimate 



\\D t U(t,-)\\ q + \\^A(D)U(t,-)\\ q + \\e(t,-)\\ q 

< (1 + t)-m-^ (||E/i|U + i + \\U 2 \\ p , r + INU) (1-7) 



for dual indices p € (1, 2], pq = p + q, and Sobolev regularity r > 2(l/p — l/q). 

If the coupling functions vanish to higher order we have to relate their vanishing order I to 
the order of contact 7 between the Fresnel curve and its tangent in the corresponding direction 
and for the corresponding sheet. In Theorems 13.51 and 13.61 we show that in this case the 1/2 in 
the exponent is changed to 1/ mm(2£,7). 

Our main motivation to write this paper is to provide a unified way to treat anisotropic 
models of thermo-elasiticity. New analytical tools presented in these notes generalise to higher 
dimensions and allow to treat especially models in 3D (outside degenerate directions). The two- 
dimensional results from jj] for cubic media and [1] for rhombic media are contained / extended; 
general anisotropic media can be treated. This is discussed in some detail in the second part |16| 
of this note. 

Acknowledgements. The first author thanks Prof. Wang Ya-Guang (Shanghai Jiao Tong 
University) for the discussion about some basic ideas of the approach presented in this paper 
during his stay at the TU Bergakademie Freiberg in August 2004. The stay was supported by 
the German-Chinese research project 446 CHV 113/170/0-2. 



We use a partial Fourier transform with respect to the spatial variables to reduce the Cauchy 
problem for fUTJ) to the system of ordinary differential equations 



parameterised by the frequency variable £. 

We denote by yc\^X2 £ C 00 ^ 1 ) the eigenvalues of A(rf) and by ri,r2 € C 00 (§ 1 ,§ 1 ) corre- 
sponding normalised eigenvectors. Both depend in a real-analytic way on rj G S . In a first step 
we reduce (|2.ip to a first order system. For this we use the diagonaliscr of the clastic operator, 
i.e. the matrix M(n) ~ (ri (77)^2(77)) build up from the normalised eigenvectors, and define 



2. General treatment of the thermo-elastic system 



J7 4t + |C| 2 A(77)J7 + i7^ = 0, 

e t + ^\ 2 e + i^-u t = 0, 

U(0,-) = U u U t (0,-) = U 2 , 0(O,-) = 0o 



(2.1a) 
(2.1b) 
(2.1c) 



UW(t,0=M T (r))U(t,0. 



(2.2) 
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Then we define by the aid of 

v 1 / 2 ( v ) = di ag (cj 1 (v),Mv)), ^(^y^erts 1 ), (2.3) 

the vector-valued function 

f(D t + v^(0)u {0) (t,0\ 

V {0 Ht,0 = (D t -VV 2 (0)UW(t,0 , (2.4) 

V J 

where as usual D t = — idf. It satisfies a first order system with apparently simple structure. A 
short calculation yields D t V(°\t,£,) = £(£)^ (0) with 

/ MO i 7 ai(6\ 

B(0= -^(0 hai(0 , (2.5) 

-w 2 (£) 1702(C) 
\-£ai(0 -|«2(0 -^oa(0 i^lCl 2 / 

where we used the coupling functions 

%(0 = ^w-e (2.6) 

For later use we introduce the notation -Bi(£) and -82(C) f° r the homogeneous components of 
-B(£) of order 1 and 2, respectively. The coupling functions aj(rj) can be understood as the 
co-ordinates of 77 with respect to the orthonormal eigenvector basis {^1(77), r 2 (j])}. Therefore, it 
holds a\(rj) +03(77) = 1. Furthermore, they are well-defined and real-analytic functions on S 1 . 

In the following proposition we collect some information on the characteristic polynomial of 
the matrix £?(£). 

Proposition 2.1. (1) tr B(£) = i/i|£| 2 and det B(£) = ik\£\ 6 dct A{rj). 

(2) The characteristic polynomial of B(£) is given by 

det(i/I - B(0) =(" - i^| 2 )(^ 2 - M0)(" 2 ~ MO) 

- vi 2 \eal(v)(v 2 ~ MO) ~ vi 2 \Z\ 2 4{v)(v 2 MO)- (2-7) 

(3) ^4n eigenvalue v £ speci?(£), £ 7^ 0, is reaZ «/ and on/y */ ^ 2 = >%o(£) / or an index 
jo = 1,2. J/ i/ie direction is non- degenerate this is equivalent to a,j (rf) = 0. 

(4) Ifaj(rj) 0, j = 1,2 £/ie eigenvalues v £ spcci?(£) satisfy 



HCI 2 , 2 a 2 (Q , 2 al(Q 

MO » 2 -MO' 



1 = + 7 2 -^L^ + 7 2 2 " 2V % fV (2.8) 



It turns out that the property of -B(£) to have real eigenvalues depends only on the direction 
77 = € S 1 . We will introduce a notation. 

Definition 1. We call a direction 77 £ S 1 hyperbolic if -B(f) has a real eigenvalue and parabolic 
if all eigenvalues of -B(£) have non-zero imaginary part. 

In hyperbolic directions we always have a pair of real eigenvalues. If 77 £ S 1 is hyperbolic 
with ±Wj (£) € spec-B(^) for £ = |£|r7, we call 77 hyperbolic lozf/i respect to the index jo (or 
with respect to the eigenvalue Xj (rj) of ^(77)) and v±(0 — i^ioCO the corresponding pair of 
hyperbolic eigenvalues of £?(£). 

A non-degenerate direction is parabolic if and only if 07(77) 7^ 0, j = 1,2, while for non- 
degenerate hyperbolic directions one of the coupling functions aj (77) vanishes. Degenerate di- 
rections are always hyperbolic (in 2D), see (|2.7p . 
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Example 2.1. If the medium is isotropic, A(rf) = /il + (A + 77)77® 77, the eigenvalues of A are fx and 
A + /i with corresponding eigenvectors 77 and r/ 1 - . Thus all directions are hyperbolic (with respect 
to the second eigenvalue). In this case the matrix decomposes into a diagonal hyperbolic 
2 x 2-block and a parabolic 3 x 3-block. This decomposition coincides with the Hclmholtz 
decomposition as used in the standard treatment of isotropic thermo-elasticity. 

Example 2.2. For cubic media (where we assume in addition 77 7^ r and /i + A 7^ 0) there exist 
eight hyperbolic directions determined by 771772 = or 77 J =772- The functions 07(77) have simple 
zeros at these directions. 

Example 2.3. Weakly coupled cubic media with A + 77 = 0, /j 7^ r, have the degenerate directions 
rjl = r/2, media with 77 = r, A + 77 7^ 0, for 771772 = 0. In both cases the coupling functions 0^(77) 
do not vanish in these directions. 

If ji = t = — A, the elastic system decouples directly into two wave equations with propagation 
speed y/Ji. In this case all directions are degenerate. 

Example 2.4. For rhombic media we have to distinguish between three cases. 

Case 1. If the material constants satisfy (A + 2/j — n)(A + 2/j — r 2 ) > 0, we are close to the cubic 

case and there exist eight hyperbolic directions given by 771772 =0 and 

77?(A + 2/j - n) = nl{\ + 2/j - t 2 ). 

Case 2. If we assume on the contrary that (A + 2/j — ti)(A + 2/i — t 2 ) < 0, only the four hyperbolic 
directions 771772 = exist. In the Cases 1 and 2 in each hyperbolic direction one of the coupling 
functions 0^ (77) vanishes to first order. 

Case 3. In the borderline case n = A + 2/j or r 2 = A + 2/j, but n 7^ r 2 , three hyperbolic directions 
collapse to one. We have the four hyperbolic directions 771772 = 0, at two of them (77^ = ±1 if 
Tj = X + 2/i) the vanishing order of the coupling function is three. 

Rhombic media are degenerate if a) 77 — T\ (or 77 = t 2 ) with degenerate direction (0, 1) T (or 
(1,0) T ) or b) A + 77 = (weakly coupled case) and (/U — Ti)(/i — T2) > with degenerate directions 
determined by 771(77 — ri) = 77K77 — T2) or c) Tj = /i = —A (exceptional case) where all directions 
are degenerate. 

Proposition 2.2. Let the direction 77 € S 1 be non- degenerate and hyperbolic with respect to the 
index jo . Then the corresponding eigenvalues v± (£) satisfy 

lim - 2 ( ,f i0 m = mm = c,t iAjiei (2.9) 

/or oZZ non-tangential limits with real constants Cfj, G R, Z?^ > 0. Furthermore, the imaginary 
part of the hyperbolic eigenvalue satisfies 

li m ^ = ^|^ >0, (2.10) 
^ a| (77) 2kC? + |£| 2 £? 

and i/ms vanishes like lm.v(\£\rf) <~ 0^(77) as n ^> f) for all £ 7^ 0. 
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Proof. Let for simplicity jo = 1. We use the characteristic polynomial of B(£) to deduce 



-+1 T ML t (2.11) 

ui(rj) xi(r?) - x 2 (r?) 

= 1- ,-, Ti^ie|-7^(iei)- (2-12) 

The existence of the limit is implied by i/± 7^ and 4 7^ as consequence of ^±(£) — > 

±|^|wi(r/) as 7/ — > 77 by Proposition ^. II 

Obviously, Im<7^(|£|) = =fk|£|/7 2 wi(?7) = =p£)^|^| is non-zero for £ 7^ and considering the 
imaginary part of the first limit expression 

a?(0 



Ini(fa(|£|) = limlm 
= lirn 



-2 Re MO Imi/±($K(Q 

l4(0-*i(0l a 



= T2^i(|£b7)k I1 (lel)| 2 lim 

17-^) |£raf(77) 

proves the second statement, lim,,^ ^ft^p = 2o Jl (»j)(c]+|ci 2 -D^) ' D 

In the case of isolated degenerate directions (others are not of interest, because then the 
system is decoupled) we can find a replacement for Proposition 12. II 

Proposition 2.3. Let fj £ S 1 be an isolated degenerate direction, Xi(fj) — ^(fj). Then the 
corresponding hyperbolic eigenvalues f±(£) satisfy 

^ "f-^f =al(n)>0, (2.13) 
n-*v uji(Q - w 2 (£) 

and, therefore, 

fl-»flwi(0-«a(0 wi(0-«2(0 

Thus, if ai(fy) 7^ then the eigenvalues v±(£) approach ±u;i(£) at the contact order between 
wi(£) and w 2 (£) (while they approach ±u>i(£) with a higher order if ai(£) = 0). 

2.1. Asymptotic expansion of the eigenvalues as |£| — > 0. If |£| is small the first order part 
Bi(£) dominates B 2 (£)> so the properties of the eigenvalues are governed by spectral properties 
ofS x (0. 

Proposition 2.4. (1) trBi(£) = and detBi(^) = 0. 

(2) If the direction r\ G S 1 is parabolic the nonzero eigenvalues v of Bi(r/) satisfy 

v z - xi(ri) v z - x 2 [ri) 
and are thus real and related to Xj(rf) by 

< xi fa) < %(n) < Mri) < %(ti) (2.16) 

(ifxi(ri) < x 2 {rj)). 
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(3) If n is non- degenerate and hyperbolic with respect to the index 1 we have — vf(r)), 

while for hyperbolic directions with respect to the index 2 three cases occur depending on 
the size of the coupling constant 7: 

V < ^{fj) - xi(fj) : ^{rf)=vl{rj), 

I 2 = X2(fj) - xiifj) : i>l{rf) = x 2 (n) = D%(v), ( 2 - 17 ) 
j 2 > x 2 {ff) - Ki(rj) : Ox(v)= ^(v)- 



(4) If the direction is degenerate, — Kyfj]), we have the eigenvalues ± -y/ >c\ (77) and 

± V >( i(il) +7 2 - 

The existence of five distinct eigenvalues of the homogeneous principal part B\(rf) for all 
parabolic and most hyperbolic directions allows us to calculate the full asymptotic expansion of 
the eigenvalues v(£) of B{£) as |£| — > 0. We will give only the first terms in detail, but provide 
the whole diagonalisation procedure. Assumption (A4) guarantees the non-degeneracy of #i(£). 
Note that, even if (A4) is violated the matrix is diagonalisable (as consequence of its block 

structure) . 

Proposition 2.5. As |£| — > the eigenvalues of the matrix -B(£) behave as 

M0=^\X(v)+O(\(i\ 3 ), (2.18a) 
i/f (0 = ±m(v) + i*\S\ 2 bM + 0(\tf) (2.18b) 

for all non-j- degenerate directions, where the functions bj € C°°(S 1 ) are given by 

W=(l + ^¥ + ^f¥)">0 (2.19) 



and 



t'A'D- h - v J "f''/' .-'i . 1 >o. (2.20) 



Furthermore, bj (77) > if r] is parabolic and bj (rj) — Q if n is hyperbolic with respect to the index 
Jo- 

Proof. We apply a diagonalisation scheme in order to extract the spectral information for £?(£). 
We assume that the eigenvalues are denoted such that *c\{rj) < x^iv)- 



Step 1. By Proposition 12.41 we know that the homogeneous first order part Bi(rj) has the 
distinct eigenvalues Dq — and ^fiv) = ^jiv)^ which are ordered as xi{n) < v\{rj) < x%{y\) < 
V2 (77) (where equality holds only under the exceptions stated in Proposition 12. 4p . We denote 
corresponding normalised and bi-orthogonal left and right eigenvectors of the matrix Biirf) by 
je ± (?7) and e^{rj). If we collect them in the matrices 

L(V) = (oe(v)\ie + (v)\ie-(v)\2e + (v)\2e-( V )), (2.21a) 

R(V) = (e Q (vM(v)\^(vM(vM(v)), (2.21b) 
we have L*[rj)R{rj) = I and 

L%r J )B 1 (r,)R( V ))=V 1 (ri)=dia,g(0,h(v),-h(v),Mv),-h(v))- (2.22) 

Further we get 

L* (v)B 2 (v)R(v) = ™h (77) ® jffi, (2.23) 
where 6* (77) and *b(n) are vectors collecting the last entries bo(rf), b^{rj) and ob{rj), ^(rf) of the 
eigenvectors 60(77), ef(n) and o^iv): j e± (v)i respectively. 
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The matrix 

B<°>(0 = L*( V )B(OR(v) = \Z\Vi(v) + \e^( V ) ® Z&fa) (2.24) 
is diagonalised modulo 0(|£| 2 ) as |£| — ► and has a main part with distinct entries. We denote 
i?( 2 )(£) = \tfi K K®7b. 

Step 2. We construct a diagonaliser of B^(^) as |£| — > of the form 

iVfe(0=^ + El^' iVO ' ) ('?)- (2-25) 

For this we denote the fc-homogeneous part of i?( fc )(£) by R( k \£) and its entries by R.[j\r]). 
Then we set for fc = 1, 2, . . . 

P fe+1 (7?)=diag J R( fc+1 )(r / ), (2.26) 



d 2 (r))-di(r;) d 2 (r))-d5(r;) 



\d 5 (r))-di(r,) d 5 ( V )-d 2 {v) U / 



(2.27) 



where g?j (77) are the entries of T>i(rj). By construction we have the commutator relation 

[V 1 ( V ),N( k \r ] )}=V k+1 (r 1 )-M k + 1 \ V ), (2.28) 

such that 

fe+i 

^+2)^ = S (°)(o^ fe (o - iv fe (o e i^i^-w 

= i? (fc+D (?) 

fc+1 

+ iei fe s (0) (^ (fc) w - ici fc ^ (fc) w E \^ v m N k m\ k+l v k+1 ( V ) 

i=i 

fe+i 

= o(iei fe+2 ). 

Using Nk(£,) — I = 0(\£\) we see that for |£| < Ck, c k sufficiently small, the matrix Nk(£) is 
invertible and 

fc+i 

N^iOBMWM) = E + G(|e| fc+2 ). (2.29) 

Thus, the entries of / Dj(rj) contain the asymptotic expansion of the eigenvalues, while the rows 
°f Nk(£,)R(r]) (and L^NjT (£)) give asymptotic expansions of the right (and left) eigenvectors 
ofS(0- 

Furthermore, the construction implies that all occurring matrices are smooth functions of 

r\ e S 1 . 

5£ep 5. We calculate the first terms explicitly. For this we need the diagonal entries of the 
matrix <g> *6. Therefore, we determine the left and right eigenvectors of B\{rj). If we assume 
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that the direction 77 is non-degenerate we get for eo(r)) = (r 1 f ,r 1 ,r% ,r 2 , r ) T and o e (??) = 
(£f , £1 , £2 , £2 , £o) T the equations 

±rfujj + vyajro = 0, ±ifuj - ^aX = 0, (2.30a) 

2 

airf + a 2 r^ + airf + a 2 ^ = 0, a,\£\ + a 2 £f + + #2^2 = 0; (2.30b) 
together with the normalisation condition 

rt~tt + ---r^ + r T Q = 1. (2.31) 
The first equations imply the representation 



J J 2^(77) 

the second line of equations follows from the first, while the normalisation condition yields 

To calculate the eigenvectors we can further require ro(r)) = £o(rj) = \fbo{rj) > 0. 

Similarly, we obtain for the eigenvectors e^" (77) = (rf , . . . , , ro) T and fce + (?7) = . . . , £2, £o) T 
the equations (we use the same notation as above in the hope that this will not lead to confusion 
here) 

±rfuj 3 + i iaj r = ±v k rf , ±£~f^ - ^4 = ±h (vjtf (2-34) 

together with the normalisation condition. Thus for parabolic directions we get 



±rf(v)= ~ , \ J ( , Mv), ±£f(r!) = - ' JK , (2-35) 



and hence 



62 fo) - r („) w = (1 + E 2 + E 2 2^T^w J 

For e~^(jf) and ke-~{rf) we have to replace 1^(77) by —Vkiv) and obtain = 6^ (??) = bk(i]). 

If the direction 77 is non-degenerate and hyperbolic with respect to the index jo, the entries rf o 
and £^ Q are undetermined by (|2.34p . while the other entries of the vectors are zero. Together with 
the normalisation condition this determines the eigenvectors and gives bj (rj) — 0. It remains to 
consider degenerate directions. Then we have v\ = oj\ and i>2 > io\ such that for k = 1 we have 
£0 = I'd = 0, rf and if are non-zero while = if — 0, especially 61(77) = 0. For k — 2 we get 
from the above expression for 62(77) = (2 + 2xr 2 /7 2 )~ 1 > 0. □ 

Remark. 1. Note that for all non-degenerate hyperbolic directions 77 6 S 1 with respect to the 
index 1 the limit 

lim ^(77)6^(77) = -1*107) (l - —f^ 1 j-c) (2-37) 

is taken and non-zero, while for hyperbolic directions with respect to the index 2 the corre- 
sponding limit is non-zero only if 7 2 7^ — Ki(fj), i.e. if the direction is not 7-degenerate. 
Near 7-degenerate directions Step 1 of the previous proof is still valid. Similar to Step 2 we 
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can diagonalise to a (2, 2, 1) block structure. The eigenvalues of these blocks can be calculated 
explicitely. 

2. For degenerate directions we obtain similarly 

lim = a i^) Q 2^) (2.38) 

v->fi (xi(r/) - X2{ti)) 2 2^ 2 xi{f]) 

and 61(77) vanishes to the double contact order. 

2.2. Asymptotic expansion of the eigenvalues as |£| — ^ 00. If we consider large frequencies 
the second order part -62(C) dominates i?i(£). This makes it necessary to apply a different two- 
step diagonalisation scheme. We follow partly ideas from [9], [H], [15] adapted to our special 
situation. 

Proposition 2.6. As |£| — + 00 the eigenvalues of the matrix B(£) behave as 

M0=te\tf-— + O(\S\~ 1 ), (2.39a) 

K 

if (£) = ±\t\u>M + %a*{r,) + Oder 1 ). (2.39b) 
for all non- degenerate directions £/|£| € S 1 . 

Remark. 1. Note that, while in hyperbolic directions we always have *f (£) = ±|£|wj (»7) for 
one index jo, in degenerate hyperbolic directions all 07(77) may be non-zero. Hence the statement 
of the above theorem cannot be valid in such directions in general. 

2. Degenerate directions play for large frequencies a similar role as 7-degenerate directions play 
for small frequencies. 

Proof. The proof will be decomposed into several steps. In a first step we use the main part 
B 2 (0 = i^|£| 2 diag(0, 0, 0, 0, 1) to block-diagonalise B(£). In a second step we diagonalise the 
upper 4x4 block for all non-degenerate directions. 

Step 1. For a matrix B € C 5x5 we denote by b-diag 4 x B the block diagonal of B consisting 
of the upper 4x4 block and the lower corner entry. We construct a diagonalisation scheme to 
block-diagonalise £?(£) as |£| — > 00. 

We set i? (_1) (£) = B(£) - B 2 (0 = and B_ 2 (0 = B z{0 and construct recursively a 

diagonaliser modulo the upper 4x4 block, 

k 

M k (Z)=I + J2\£\~ j MW(7 1 )- (2-40) 
i=i 

Again we denote by Ry^(rf) the (— /c) -homogeneous part of (which exists because it exists 

for (£) and the existence is transfered by the construction). Then we introduce the recursive 
scheme 

B k - 2 ( V ) =b-diag 4ilJ R( fc - 2 >(r,), (2.41) 



mW(ti) = - 







\-R ( t 2) (v) ■•• -R%T 2) (v) 

for k = 1, 2, . . ., such that the commutator relation 

[B-aM.AfWfo)] =B fc _ 2 (7 ? )-^- 2 )(r ? ) ) (2.43) 



^ 2) (^7) 



(2.42) 
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holds. Thus it follows 



(2.44) 



fc-2 

R< k - 1 HO=m)M k (t)-M k ® £ \t\- j B j (r 1 ) = 0(\e- k ) 

and using that Mfc(£) is invertible for |£| > Cfc, Ck sufficiently large, we obtain the block diago- 
nalisation 

fc-2 

M-\oBmuo = Yl \^r^M+o(\e- k ), (2.45) 

where Bj(r/) = b-diag 41 /Sjfa) is (4, l)-block diagonal. 

Step jg. By Step 1 we constructed M k (Q such that x (£).B(£)Mfc(£) is (4, l)-block diagonal 

modulo 0(|£| 1_fc ). The upper 4 x 4 block has already diagonal main part | ^ | 1 X? 1(^7) = -Bi(£)- ^ 

the direction 77 = £/|£| is non-degenerate, the diagonal entries iwifa) and iu>2(v) are mutually 
distinct and thus we can apply the standard diagonalisation procedure (cf. proof of Proposi- 
tion 12. 5p in the corresponding subspace. This does not alter the lower corner entry and gives 
only combinations of the entries of the last column and of the last row (without changing their 
asymptotics). 

Thus we can construct a matrix N k -x(g) = I + T,jZl N^^), which is invertible for 
|£| > Cfc_i, Cfe_i sufficiently large, such that 

fc-2 

N^iOM^iOBiOMkiONk-iiO = \eB-2(v) + E Itl-'VM+Om 1 -") (2-46) 

i=-i 

is diagonal modulo 0(|£| 1-fc ). 

Step 3. We give the first matrices explicitly Following Step f we get 

10-2(11) 

7^2(17) 



-yaiiv) 7 Q 2(i7) r yai(v) 7 Q 2(i)) 
2k 2k 2k 2k 



(2.47) 



/ 



together with B-\{rf) = diag(wi(7j), o^fa), — ^1(77), — u>2(t])i 0) an d 



( jT^M 



.7 a 1 (ri)a 2 {ri) 



2k 

-l' A ai(ri)a 2 (v) 
2k 



7 2 affa) 
2k 



701(17)^1(1)) 
2k 



7a 2 (T))^2(')) 
2k 



701(17)^1 (17) \ 
2k > 

702(17)^2(1)) 
2k 



i2i y 



(2.48) 



Bo fa) =b-diag 4 . 1J R(°)(7 ? ) 



(2.49) 

in the first diagonalisation step. Applying a second step alters only the last row and column to 

OdCr 1 ). Following Step 2 we diagonalise the upper 4x4 block to |£|S-i(r?) + B (rf) 

modulo 0(|£| _1 ) and the statement is proven. □ 

Remark. If the direction 77 is degenerate, i.e. ^1(77) = Wafa), we can block-diagonalise in Step 
2 to (2, 2, 1)-block form. To diagonalise further we have to know that the 0-homogeneous part 
of these 2 x 2-blocks has distinct eigenvalues. 



12 



MICHAEL REISSIG AND JENS WIRTH 



One possible treatment of degenerate directions is given in the following proposition. Note 
that a corresponding statement can be obtained for 7-degenerate directions as |£| — > 0. 

Proposition 2.7. Let fj be an isolated degenerate direction. Then the corresponding hyperbolic 
eigenvalue satisfies in a small conical neighbourhood of fj 

jo ^ } 2 4k 



MiH^iZ ^ 1 i7 2 M0-"2(Q)(a?fo)-a!fo)) 
4 16k 2 4k 

+ 0(|er 1 ). (2.50) 

Proof. We follow the treatment of the previous proof to (2, 2, l)-block-diagonalise -B(£) modulo 
Now, we consider one of its 2 x 2-blocks. (We use a similar notation as before in the hope 
that it will not lead to any confusion.) Such a block is given by 

8(0 = IflB-ifa) + 80(7?) + Oder 1 ), (2.51) 

where 

B_i(i 7 ) = diag(wi(t ? ),w 2 (»7)) 1 (2.52) 

Bbfo) = !£( , ^(^)V (2 . 53) 

In the direction 77 both diagonal entries of B-% coincide. In a small conical neighbourhood we 
denote the eigenvalues of |£|£?-i(?7) + Bo{rj) as 6+(£) and £-(£)• A simple calculation yields 

, , c , ^i(0+^2(0 , i7 2 

M0 = 2 + ^ 



± V 4 16^ + 4^ (2 ' 54) 

with <5_(£) = oji{C) and 5+ (£) = + fer. The hyperbolic eigenvalue corresponds to £_(£)• 

These eigenvalues are distinct in a sufficiently small neighbourhood of fj (and may coincide only if 
01(77) = a 2 (?7) and —i^ziO) 2 ~ 7 4 /(4 k2 ) ; which gives eventually two parabolic directions). 

Hence the perturbation theory of matrices implies = + C ( | <^ | 1 ) in a sufficiently 

small neighbourhood of 77 and the statement is proven. □ 

Remark. Note, that 

for all fixed which coincides with the result (|2.13|) of Proposition 12. 31 for the eigenvalue 



2.3. Collecting the results. The asymptotic expansions from Propositions 12.51 and Proposi- 
tion 12.61 imply estimates for eigenvalues of B(£) and the proofs give representations of corre- 
sponding eigenvectors. For the application of multiplier theorems and the proof of L p -L q decay 
estimates it is essential to provide also estimates for derivatives of them. 

Assume that the eigenvalues under consideration are simple. From the asymptotic expansions 
we know that this is the case for small frequencies and also for large frequencies. For the middle 
part it will be sufficient to know that the hyperbolic eigenvalues are separated, which follows for 
sufficiently small conical neighbourhoods of these directions. 
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In a first step we consider derivatives of the eigenvalues. Differentiating the characteristic 
polynomial 

5 

= det(*/(0/ - B(0) = E WHO" (2-56) 

k=0 



with respect to £ yields by Leibniz formula 

» = EE& D ^K))( D (^) (2-57) 



for all multi-indices a <E Ng. Thus we can express the highest derivative of v(S,) k in terms of 
lower ones and hence Faa di Bruno's formula (see e.g. [5]) yields an expression 

(dmo) E^KMO* -1 = E E c k ^(n a -?i k (0)(D^(0) (2-58) 

fe=l fc=0/3<a 

with certain constants Ck,a,/3- Because the eigenvalue has multiplicity one, the sum on the left- 
hand side is nonzero and therefore we can calculate the derivatives of z/(£) by this expression. 
Furthermore, it follows that for small and large frequencies the derivatives of the eigenvalue have 
full asymptotic expansions and thus we are allowed to differentiate the asymptotic expansions 
term by term. 

It remains to consider the corresponding eigenprojections. Recall that if is a eigenvalue of 
multiplicity one and r(£) and are corresponding right and left eigenvectors, the correspond- 
ing eigenprojection is given by the dyadic product P^iO = l{0 ® r (£)- Thus the constructed 
diagonaliser matrices imply asymptotic expansions of these operators. Again we are only inter- 
ested whether the derivatives of these eigenprojections also possess asymptotic expansions (in 
order to see whether it is allowed to differentiate term by term). 

For this we use the representation 

^(o= n WQi-Bmm-v®)- 1 (2-59) 

j>especB(5)\{f} 

given e.g. in [I], [7]. All terms on the right-hand side have full asymptotic expansions as |£| — > 
and |£| — > oo together with all of their derivatives. Differentiating with respect to £ yields the 
same result for the eigenprojection. Thus we obtain 



Proposition 2.8. The asymptotic expansions from Provosition \2.5\ and Provosition \2. 61 may be 
differentiated term by term to get asymptotic expansions for the derivatives of the eigenvalues. 
Furthermore, the same holds true for the corresponding eigenprojections. 

From Proposition 12 . II we know that an eigenvalue v(t;) of the matrix B(£) is real if and only if 
rj = £/|£| is hyperbolic. We want to combine this information with the asymptotic expansions of 
Proposition ^. 51 and Proposition ^. 61 and derive some estimates for the behaviour of the imaginary 
part. 

Proposition 2.9. Let c > be a given constant. 

(1) Let n — £/|£| € S 1 be parabolic. Then the eigenvalues of B(£) satisfy 

hnvfe) > C > 0, |f|>c, (2.60) 
with a constant C n depending on the direction r\ and c. Furthermore, 

lmv(0~ b (r))\e, |£|<c, (2.61) 
where b(rj) is one of the functions from Proposition[ 
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(2) Let fj be non-degenerate and hyperbolic with respect to the index 1. Then (\t;\fj) = 
±\£\u>i(fj) and Vo(£) o,nd u 2 (£) satisfy the statement of point 1. Furthermore, 

luivf{t)~a 2 M, \£\>c, \ri-fj\^l, (2.62) 

and 

lmvf(0~\eal(v), lel<c, \V-V\<Z1. (2.63) 

Proof. The first point follows directly from the asymptotic expansions, we concentrate on the 
second one. We know that the hyperbolic eigenvalues vf (£) satisfy by Proposition 12.21 

Imvt(0 = 4(v)Nf(0 (2.64) 

with a smooth non- vanishing function (£) defined in a neighbourhood of fj. By Proposition ^. 51 
and 12.61 we see that N^^) also has full asymptotic expansions and thus 

Nf t (0 = ^+O(ICr 1 ), |£|-oo (2.65a) 
Nf (0 = H£| 2 |j| + ^Cl€| 3 ), ICI - 0. (2.65b) 

Together with (|2.37p we get upper and lower bounds on N^(^) and the desired statement follows. 

□ 

Remark. A similar reasoning allows to replace the hyperbolic eigenvalue near degenerate di- 
rections by the model expression obtained in Proposition 12.71 thus ~ uniformly in a 
sufficiently small conical neighbourhood of fj. 

3. Decay estimates for solutions 

Our strategy to give decay estimates for solutions to the thermo-elastic system (|l.ip is to 
micro-localise them. In principle we have to distinguish four different cases. On the one hand 
we differentiate between small and large frequencies, on the other hand between hyperbolic 
directions and parabolic ones. 

We distinguish between two cases depending on the vanishing order of the coupling functions. 
If the coupling functions vanish to first order at hyperbolic directions only, we rely on simple 
multiplier estimates. Later on we discuss coupling functions with higher vanishing order, where 
the decay rates are obtained by tools closely related to the treatment of the elasticity equation. 

3.1. Coupling functions with simple zeros. In a first step we consider the first order system 

B t V = B(D)V, V(0,-) = V (3.1) 

to Cauchy data V <E 5(R 2 ,C 5 ). For a cut-off function x € C°°(M + ) with x(s) = 0, s < e, and 
x(s) = 1, s > 2e, we consider 

P P ar(v)= II X(\V-fj\), Phyp{ri) = l-P P ar(ri). (3.2) 

fj hyperbolic 

Then Ph yp (D) localises in a conical neighbourhood of the set of hyperbolic directions, while 
Ppar(D) localises to a compact set of parabolic directions. 

The asymptotic formulae and representations for the characteristic roots of the full symbol 
.£?(£) allow us to proof decay estimates for the solutions. 
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Theorem 3.1. Assume that (Al) to (A4) are satisfied and the coupling functions vanish in 
hyperbolic directions to first order. 

Then the solution V(t, x) to (|3.1|) satisfies the following a-priori estimates: 

\\ X (V)Ppar(V)V(t, -)\\ q < e- Ct ||Vo|| Pir (3.3a) 

||(1 - x(D))P por (D)F(i, .)||, < (l +t)-G-$\\V \\ p (3.3b) 

||x(D)i\ OT ,( D )^(*.-)H a ^(l + *)" i( *"« ) l|V r o||p,r (3.3c) 

||(1 - X (D))P fcw (D)F(i, OH, < (1 + t)-W-V\\V \\ p (3.3d) 
for dual indices p S (1, 2], pq = p + q, and with Sobolev regularity r > 2(l/p — 1/q). 

Remark. If B(£) is diagonalisable for £, ^ (which is valid e.g. if P(f ) has no double eigenvalues 
for £ ^ 0) we can make the result even more precise. To each eigenvalue G specP(£) we 
have corresponding left and right eigenvectors and associated to them the eigenprojection Pj,(D) 
such that 

P,(D)y(i,-) = e iMD) P,(D)y . (3.4) 

Thus, we can single out the influence of one eigenvalue in this way. Note, that this is only of 
interest in the neighbourhood of hyperbolic directions and for the corresponding eigenvalue and 
there the assumption of diagonalisability of P(£) may be skipped (real eigenvalues are always 
simple, for small |£| diagonalisation works, for large |£| everything goes well under the assumption 
of non-degeneracy and on the middle part we make the neighbourhood small enough to exclude 
possible multiplicities). 

Proof. We decompose the proof into four parts corresponding to the four estimates. The micro- 
localised estimates are merely standard multiplier estimates. We do not use stationary phase 
method. 

Step 1. Parabolic directions, large frequencies. In this case we have uniformly in £ S supp xPpar, 
the estimate Im z/(£) > C > 0. Taking < C < C we obtain for these £ 

*nB(Z) = B ®- 2 . B '®>CI 
in the sense of self-adjoint operators and the estimate follows from the L 2 -L 2 estimate 

^\\x(P)Ppar(Pmt,x)\\l = ^ t \\x(0P P ar(0V(t,0\\ 2 2 
= 2 Re (x(0P P ar(0V(t, 0, X(0P P ar(mV(t, £)) 

= -2 Im (x(Z)P P ar(0V(t, 0, x(0P P aA0B(0V(t, 0) 

< -2C\\ X (0Ppar(t)V(t,0\\ 2 = -2C\\ X (D)Ppar(P)V{t,x)\\l 

viewed as H s -H s estimate and combined with Sobolev embedding. 



Step 2. Parabolic directions, small frequencies. We know from Proposition 12 . 51 that in this case 
the matrix has only simple eigenvalues. We will make use of the representation of solutions 

as 

V(t,x) = J2 j tv(V) Pu(V)V (x) (3.5) 

^)espccB(C) 

with corresponding eigenprojections (amplitudes) P u (0- The amplitudes are uniformly bounded 
on the set of all occurring £ and possess full asymptotic expansions in |£| as £ — * together with 
their derivatives. Especially, by Hormander-Mikhlin multiplier theorem 10, p. 96, Theorem 3] 
the operators P^(D) are L p -bounded for 1 < p < oo. 
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It remains to consider the model multiplier e 1 *"^'. From Proposition 12 . 51 we know that 

\tP*<&\<e- c W, |f|<c, (3.6) 
with suitable constants c and C and thus the L 1 -^ 00 estimate 

I| e itl/(D) /Iloo < He^/lli < \\^\\m m <c})\\Hco 

<n/iii fe-^ieidiei^a+r 1 ii/iu 

Jo 

holds for all / e I-^R 2 ) with supp/ C {|£| < c}. Riesz-Thorin interpolation [TJ Chapter 4.2] 
with the obvious L 2 -L 2 estimate gives the desired decay result. 

Step 3. Hyperbolic directions, large frequencies. We take the conical neighbourhoods small 
enough to exclude all multiplicities (related to the eigenvalue which becomes real in the hyperbolic 
direction). Then similar to (13. 5|) the solution is represented as 

V(t,x) = e it ^P+(B)Vo(x)+e- it ^p-(B)V (x) + V(t,x), (3.7) 

where V(t,x) corresponds to the remaining parabolic eigenvalues of and satisfies the es- 
timate from Step 1. Again we can use smoothness of P^{£,) together with the existence of a 
full asymptotic expansion as |£| — > oo to get L p -boundedness of P^(D) for 1 < p < oo from 
Hormander-Mikhlin multiplier theorem. 

It remains to understand the model multiplier e ±ltv ^ for the hyperbolic eigenvalue related 
to the hyperbolic direction fj. Using the estimate from Proposition 12.91 we conclude for r > 2 

lle^^/lloo < We^^fh < \\e ±H ^\- r \\ LHSl )\\ ItffU 

/oo re 
\e~ r m j e- c ^ 2 d<t> 

<t- 1/2 ||<D) V|| i, t>\ 

for all / S (D) _r i 1 (R 2 ) with supp/ C Si = {|£| > c, \<q — fj\ < e}. Riesz-Thorin interpolation 
with the L 2 -L 2 estimate gives the desired decay result. 

Step 4- Hyperbolic directions, small frequencies. Like for large hyperbolic frequencies we make use 
of the representation (|3.7p to separate hyperbolic and parabolic influences. As in the previous 
cases the existence of full asymptotic expansions imply that the projections Pj t (D) are L p - 
bounded for 1 < p < oo. 

It remains to understand the model multiplier e ±xtv ^ . Using Proposition ^. 9l we have \ e ±ltu ^ \ < 
e —Cit\i\ 4> such that after introducing polar co-ordinates 

||e ±a "< D >/l|oo < We^^fh < \\e ±ttv ^\\^ Sa )\\fh 

~ 11/1,1 [ £ e ~ C2t * 2m2d ^m 

^t-^H/lli < *~ 1/a H/IU, t>l 

Jo 

for all / € L X (R 2 ) with supp / C 52 = {|£| < c, \rj — fj\ < e}. Riesz-Thorin interpolation with 
the L 2 -L 2 estimate gives the desired decay result. □ 

Remark. 1. In non-degenerate hyperbolic directions where the coupling function vanishes to 
order £ (cf. Example 12.41 case 3) we obtain by the same reasoning the weaker L p -L q decay rate 

\\x(P)Ph W (P)V(t r )\\ q <(l + t)-M$-y\\V \\ p ,r. (3.8) 
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It remains to understand whether this weaker decay rate is also sharp or whether an application 
of stationary phase method may be used to improve this. We will discuss this in Section [3.21 
2. We can extend the estimate of Theorem 13.11 to the limit case p = 1, if we include the eigen- 
projections P v (£) into the considered model multiplier and use just their boundedness instead 
of Hormander-Mikhlin multiplier theorem. This is possible, because all the multiplier estimates 
were based on Holder inequalities. 

So far we understood properties of solutions to the transformed problem (|3.1[) for the vector- 
valued function V(t, x) given by (|2.4p as 

/(D t +D 1 /2( D ))C/(°)(t ) x)\ 
V(t,x) = (Dt-V^mu^^x) , (3.9) 

V / 

where EA°)(t, x) — M T (D)U(t, x) is the elastic displacement after transformation with the diag- 
onaliser of the elastic operator. Because M(rj) is unitary and homogeneous of degree zero this 
diagonaliser is L p -bounded for 1 < p < oo with bounded inverse. Thus we have 

D t U(t,x)=M(D)(j ? I ? tyv(t,x) (3.10) 

y/A^)U(t,x)=M(D)Q i ? -3 \ tyv(t,x) (3.11) 
such that as corollary of Theorem 13.11 we obtain 

Corollary 3.2. Assume (Al) to (A4) and that the coupling functions vanish of first order in all 
hyperbolic directions. Then the solution U(t, x) and 6{t, x) to (|1.1[) satisfy the a-priori estimates 



\\B t U(t,-)\\ q + \\^A(D)U(t,-)\\ q + \\0(t,-)\\ q 

< (1 + f)"i(?-i) {\\Ux\\ p , r+1 + \\U 2 \\ P , r + \\0o\\ P ,r) (3-12) 
for dual indices p S (1, 2], pq = p + q, and Sobolev regularity r > 2(l/p — 1/q). 

Remark. Including the diagonaliser M(£) into the model multipliers of Theorem 13.11 allows to 
overcome the restriction on p and to extend this statement up to p — 1. We preferred this way 
of presenting the results because they allow to decouple both statements. Theorem 13.11 gives 
deeper insight into the asymptotic behaviour of solutions than Corollary 13 . 21 does. 

3.2. Coupling functions vanishing to higher order. We have seen that under the assump- 
tion that the coupling functions a,j (ry) S C°° (S 1 ) , dj (rj) — r) ■ rj (n) related to the symbol of the 
elastic operator A(rf) have only zeros of first order, we can deduce L v -L q decay estimates without 
relying on the method of stationary phase. 

Now we want to discuss how to use the method of stationary phase to deduce decay estimates 
in the remaining cases. From Proposition ^. 21 we know that in hyperbolic directions the imaginary 
part Im (£) of the hyperbolic eigenvalues vanishes like the square of the corresponding coupling 
function 0^(77), while the real part Re;/^(£) is essentially described by ±Wj (£). First, we make 
this more precise and formulate an estimate for Rei^ (£) =F ^joiO anc ^ ^s derivatives. 

Proposition 3.3. Let fj S S 1 be hyperbolic with respect to the index jo and a,j (r]) vanish to 
order £ in fj. Then there exists a conical neighbourhood of the direction fj such that for all 
k = 0, 1, . . . , 11 — 1 the estimates 

-l|^(Re^(0T^ (0)| <c(\mv-fj\ U - k , (3.13a) 

|fl*Imi^(0|<c(|C|)|„-^-* (3.13b) 
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hold uniformly on it, where c(|£|) ~ 1 as |£| — > oo and c(|£|) ~ |£| as |£| — > 0. 

This estimate may be used to transfer the micro-localised decay estimate for the elasticity 
equation to the thermo-elastic system on a sufficiently small conical neighbourhood of fj. We 
need one further notion to prepare the main theorem of this section. Decay rates for solutions to 
the elasticity equation depend heavily on the order of contact between the sheets of the Fresnel 
surface and its tangents, cf. :l3] or point 2 of the concluding remarks on page I2T1 

Proposition 3.4. For t] G S 1 we denote by *fj(r)) the order of contact between the j-th sheet 

= {^(77)77 1 77 eS 1 } (3.14) 

of the Fresnel curve and its tangent in the point ujj l (rj)rj. Then 

d*{d 2 v uM +w J -fa)) = > k = Q,..., 7,(77) -2, (3.15) 

(tfijiv) > V and 

tf'MuM + ^ ilv) ~ 2 uM + 0. (3.16) 

Proof. Follows by straight-forward calculation. The curvature of the j-th Fresnel curve Sj at the 
point u)~ (77)77 factorises as Wj(rj) + d 2 u>j(rj) and a smooth non- vanishing analytic function. □ 

For the following we chose the conical neighbourhood of the hyperbolic direction fj (with 
respect to the index jo) small enough in order that the curvature of Sj vanishes at most in fj. 

Theorem 3.5. Let fj be hyperbolic with respect to the index jo and let aj (77) vanish in fj to order 
£ > 1. Let us assume further that */j (fj) < 21. Then 

|| e ±i^ (D) / || 9 < (l+t) _ ^ ( *-* ) ||/llp,r (3.17) 

for dual indices p £ [1,2], pq = p + q, regularity r > 2(1/ p — 1/q) and for all f with supp/ 
contained in a sufficiently small conical neighbourhood of fj. 

Proof. We distinguish between two cases related to the Fourier support of /, for general / we 
can use linearity. In both cases we apply the method of stationary phase. Key lemma will be 
the Lemma of van der Corput, cf. [TTJ p. 334], combined for large frequencies with a dyadic 
decomposition. It suffices to consider t > 1, the estimate for t < 1 follows from the uniform 
boundedness of the Fourier multiplier together with Sobolcv embedding using the regularity 
imposed on the data, r > 2(1/ 'p— 1/q). 

In the following wc skip the index jo of the eigenvalue. Thus ^(77) stands for ujj (rf) and 7(77) 
for 7 Jo (77) to shorten the notation. 

Large and medium frequencies. Assume that / is supported in a sufficiently small conical neigh- 
bourhood of 77 bounded away from zero, |£| > 1. We will use a dyadic decomposition in the radial 
variable. Let for this \ £ Co°(BL) be chosen in such a way that xW = [0, 1], SU PPX £ [1/2,2] 
and £ 3e zX(2 J s) = 1 for all s S R+. We set Xj(s) = x(2 _J s) such that suppxj C [2 J '"\ 2^ +1 ]. 

We follow the treatment of Brenner [2] and use Besov spaces. Due to our assumptions on the 
Fourier support of / Besov norms are given by 

II/IIb;., = I|2 ,v II3o-(|d|)/IIpIL. (No) (3-i8) 

and corresponding Besov spaces are related to Sobolev and L p -spaces by the embedding relations 



- b; >2 , B% 2 ^ L q 



(3.19) 
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for p G (1, 2] and q S [2, oo). Thus it suffices to prove anBp 2 ^Bj 2 estimate for dual p and q. 
The case p = q = 2, r = is trivial by Plancherel and the uniform boundedness of the multiplier, 
we concentrate on the B\ 2 — * 2 estimate. Therefore, we use 



| e ±iMD) /l| S o 



1/2 

Ell^'( D ) e±iMD) /H-] <BUp/i(<) 



where Ij(t) are estimates of the dyadic components of the operator, 

i,{t) = sup i^i^mxMMm 



(3.20) 



(3.21) 



Here localises to a small conical neighbourhood of r\ bounded away from zero, |£| > 1, 

containing the support of /. We assume 1^(01 + l^n^iOl + l^lflVKOI ~ 1 ano - "0 = 1 on a smaller 
conic set around fj. We introduce polar co-ordinates |£| and <p (with (f> — corresponding to fj) 
and set x = tz. Using that the hyperbolic directions fj G S 1 are isolated we can assume that 
the conical neighbourhood under consideration contains no further hyperbolic direction. For the 
calculation of Ij(t) we integrate over the interval [— e, e] for 0. This yields for all j E Nq 



Ij(t) = sup 



2 3 + l 
2 3-l 



2i( 2 - r ) 



sup 

zeR 2 



1/2 



J2H(Kc v(2^)2- J +*•£) p -t Im i/(2 3 {) 



(3.22) 



Due to Proposition 13.31 (let us restrict to the — case) we have 

Reu(2^) = 2?w(() + 2 J '|£|a(2^), 
where 

\d k ^a(2^)\<C k \cf>\ 2e ~ k , k = 0,l,...,2£-l. (3.23) 

By our assumptions 2£— 1 > 7(77). We use this to reduce the problem to properties of the elastic 
eigenvalue u>j (77) and consider 



Ij(t) = 2^ 2 -^ sup 



1/2 



; i^|e|( W (r ) )+z.„+a(2^)) e - t Im,(2^) ?/ ,( 2 ^ )d ^ x( |^| ) |^|l-r de 



(3.24) 



Note that, |a(2 J £)| < \<f>\ 2e is small, if we choose the neighbourhood of fj small enough. 

If |w(?7) + 2-77! > 5 for some small S the outer integral has no stationary points and integration 
by parts implies arbitrary polynomial decay. We restrict to the z € R 2 with \u>(rj) + z ■ r)\ < 5 
and consider only the inner integral 



Ij(t,\i\) 



a it2 J |e|( W (r ) )+ 2 -T,+ Q (2^)) p -tImi/(2 3 C) 



ip{2 j £)d<f> 



(3.25) 



It can be estimated using the Lemma of van der Corput. For this we need an estimate for deriva- 
tives of the phase. We start by considering the unperturbed phase ^(77) + 2-77. Differentiation 
yields 



tfW(u>(rj) + z- V ) = d]W- 2 (dfa( V ) + u(r,)) - d^-\d^{ V ) + w(rj)) - 

\w(ri) +z-T], 



+ 



2\j(fj) 



dr,u(Tl) + Z ■ V T ', 2^(77) 



(3.26) 
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and by Proposition 13 .41 the first term is non-zero for fj while the others are small in a neighbour- 
hood of fj. Choosing 8 and the neighbourhood small enough and using (|3.23[) this implies a lower 
bound on the ^{f})-ih derivative of the phase, 

\df n) {uj(r)) + z-r) + a(2 j £))\ > 1 (3.27) 

uniformly on this neighbourhood of fj and independent of j. Thus we conclude by [111 p. 334] 
for all t > 1 and uniformly in z with \lu(t)) + z ■ r\\ < 5 

i 3 {tM\)<c k {tv\t\y^h (\ + | e |v- tIm ^ 2J ^(2 J o 

To estimate the remaining integral we use Proposition 13. 31 



(3.28) 



< 



2 r 

c 



-ct<j> 2 



< 1. 



(3.29) 

^00,2 



Integrating over £ € [1/2, 2] and choosing r > 2 we obtain for j 6 No 

;$ (i+*r^, 

where the occurring constant is independent of j and (I3.20|) implies the desired B\ 2 
estimate. Interpolation with the L 2 -L 2 estimate gives the corresponding B r v 2 — * B® 2 estimates 
with r > 2(l/p — l/q) and finally the embedding relations to Sobolev and I/ p -spaces the desired 
estimate 

l|e iMD) /ll 9 <(l+i)^ll/IU (3.30) 
for all / satisfying the required Fourier support conditions. 

Small frequencies. Assume now that / is supported in a small conical neighbourhood of fj with 
|£| < 2. In this case we estimate directly by the method of stationary phase. We sketch the main 
ideas. It is sufficient to estimate 



I(t) = sup 

iei 2 



j tv(e) mx(\(\) 



(3.31) 



where the function ip £ C°°(R+) localises to the small neighbourhood of rj with |£| < 2. Again 
we require \4>(0\ + l<V0(£)l + l^ldVKOI ~ 1- -"- n correspondence to large frequencies I(t) equals 
to 

r 8 



I{t) = sup 



„xt\£\{w{ V )+z-r)+a{£)) p -t Im u{£) 



mx{\w\i\m 



(3.32) 



We distinguish between \uj(jf) + z ■ r}\ > S, where the outer integral has no stationary points, and 
^(77) + z ■ r)\ < S. In the first case we can apply one integration by parts and get t~ l using 



(3.33) 



In the second case we can reduce the consideration to 



,it\£\{u>{ V )+z- V +a{£,)) f ,-tlmv{£,) 



(3.34) 



and an application of the Lemma of van der Corput. By Propositions 13 . 31 and I3T4I we have 



3?%fa)+*-»7 + a(O)^0 



(3.35) 
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for all £ in a sufficiently small conic neighbourhood of rj. So we obtain 

i(t, lei) ;$ weir^ (i + f |^ e - tim ^)^(e)| dA < mr^- (3.36) 

Integrating with respect to |£| yields the estimate for 7(f) 

/(<)< r/(tJCI)|£|d|£|<rw / 2 |£| 1 -w d |£|<(l + i)-^», i>l (3.37) 



./o Jo 
and by Holder inequality the L 1 -L°° estimate 

lle^/IU < /Mil/Hi < (l + r^ll/lli (3-38) 
for all / satisfying the required Fourier support condition. 

By interpolation with the obvious L 2 -L 2 estimate we get L p -L q estimates and combining 
them with the estimate of the first part proves the theorem. □ 

If the order of contact exceeds the vanishing order of the coupling funtions, the best we can do 
is to use the idea of Section f3. II and to apply standard multiplier estimates. As already remarked 
after the proof of Theorem 13.11 we obtain as decay rate in this case: 

Theorem 3.6. Let fj be hyperbolic with respect to the index jo and let aj (rj) vanish in fj of order 

1. Let us assume further that "fj (fj) > 21, where Jj (Tj) is defined in Proposition \S.4\ Then 

|| e ±it^ p) / || 9 < {1 + t) -^-i) m ^ (3 39) 

for dual indices p E [1,2], pq = p + q, regularity r > 2(l/p — 1/q) and for all f with supp/ 
contained in a sufficiently small conical neighbourhood offj. 

4. Concluding remarks 

1. In a second part of this note, [16j . we will give concrete applications of the general treatment 
presented so far. From the remarks and examples we made in the previous sections it follows that 
we indeed cover the estimates of [3] for cubic media and [1] for rhombic media in the situations 
of coupling vanishing to first order. 

In [16] we will discuss the situation of higher order tangencies for rhombic media and give a 
new estimate extending the results of [3]. Furthermore, we will give concrete examples of media 
for all achievable decay rates in the case of a differential elastic operator A(D). 

2. In general decay rates of solutions are determined by vanishing properties of the coupling 
functions. If the coupling functions are nonzero, decay rates are parabolic. If they vanish to 
sufficiently high order, decay rates are hyperbolic and determined by the elastic operator (micro- 
localised to this direction). In the intermediate case simple multiplier estimates are sufficient. 

3. In the case of isotropic media one of the coupling functions vanishes identically. In this case 
one pair of eigenvalues of the matrix -B(£) is purely real f±(£) = ±(A + /i)|£| and the components 
related to these eigenvalues solve a wave equation. Thus they satisfy the usual Strichartz type 
decay estimates, [2], 

|| e ±«(>+/*)l D lp±(D)V || g <(l + *)"* ( '~« ) ||Vo||p,r (4.1) 
for p G (1, 2], pq = p + q and r = 2(1 /p — 1/q). More generally, if we know that one pair of 
eigenvalues satisfies f±(£) = ±|£|w(?y) for all rj G S 1 with a smooth function ui : S 1 — > R+, the 
decay rates for the corresponding components depend heavily on geometric properties of the 
Fresnel curve S — { lu^ 1 ^)^ | rj € S 1 } C K 2 . Following [5] the estimate (|4.ip is valid in this case 
as long as the curvature of S never vanishes. If there exist directions rj where the curvature of S 
vanishes, the constant \ in the exponent has to be altered to 4, where 7 denotes the maximal 
order of contact of the curve S to its tangent, [T3"] . 
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4. The main focus of this paper was on the treatment of non-degenerate cases, thus assumptions 
(Al) to (A4) are required. Nevertheless, we showed how to obtain a control on the eigenvalues 
of the matrix B(£) in the exceptional cases where either (A3) or (A4) is violated. In the first 
case the treatment of large frequencies has to be replaced by the investigation of the expression 
obtained in Proposition ^. 71 while in the latter one a corresponding replacement has to be made 
for small frequencies. 

5. The results in this paper are essentially two-dimensional. For the study of anisotropic thermo- 
elasticity in higher space dimensions there arise two essential problems. The first is that we can 
not assume (A3). For example in three space dimensions and for cubic media the symbol of 
the elastic operator A : § 2 — ► C 3x3 has degenerate directions with multiple eigenvalues related 
to the crystal axes. Nevertheless, the multi-step diagonalisation scheme used in Section [521 can 
be adapted to such a case. This will be done in the sequel. The second problem is that the 
geometry of the set of hyperbolic directions becomes more complicated. In general we can not 
expect isolated hyperbolic directions, it will be necessary to consider manifolds of hyberbolic 
directions on § 2 . 
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